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ABSTRACT 

Context. The Sun traverses a low-density, hot entity called the Local Bubble. Despite its relevance to life on Earth, the conditions in 
the Local Bubble and its exact configuration are not very well known. Besides, there is some unknown interstellar substance which 
causes a host of absorption bands across the optical spectrum, called Diffuse Interstellar Bands (DIBs). 

Aims. We have started a project to chart the Local Bubble in a novel way, and learn more about the carriers of the DIBs, by using 
DIBs as tracers of diffuse gas and environmental conditions. 

Methods. We have conducted a high signal-to-noise spectroscopic survey of 670 nearby early-type stars, to map DIB absorption in and 
around the Local Bubble. The project started with a Southern hemisphere survey conducted at the European Southern Observatory’s 
New Technology Telescope and has since been extended to an all-sky survey using the Isaac Newton Telescope. 

Results. In this first paper in the series, we introduce the overall project and present the results from the Southern hemisphere survey. 
We make available a catalogue of equivalent-width measurements of the DIBs at 5780, 5797, 5850, 6196, 6203, 6270, 6283 and 6614 
A, the interstellar Na i D lines at 5890 and 5896 A, and the stellar He i line at 5876 A. We find that the 5780 A DIB is relatively strong 
throughout, as compared to the 5797 A DIB, but especially within the Local Bubble and at the interface with more neutral medium. 
The 6203 A DIB shows a similar behaviour, but with respect to the 6196 A DIB. Some nearby stars show surprisingly strong DIBs 
whereas some distant stars show very weak DIBs, indicating small-scale structure within as well as outside the Local Bubble. The 
sight-lines with non-detections trace the extent of the Local Bubble especially clearly, and show it opening out into the Halo. 
Conclusions. The Local Bubble has a wall which is in contact with hot gas and/or a harsh interstellar radiation field. That wall is 
perforated though, causing leakage of radiation and possibly hot gas. On the other hand, compact self-shielded cloudlets are present 
much closer to the Sun, probably within the Local Bubble itself. As for the carriers of the DIBs, our observations confirm the notion 
that these are large molecules, whose differences in behaviour are mainly governed by their differing resilience and/or electrical 
charge, with more subtle differences possibly related to varying excitation. 

Key words. ISM: atoms - ISM: bubbles - ISM: lines and bands - ISM: molecules - ISM: structure - local interstellar matter 


1. Introduction 

The Solar System moves through a region of unusually low gas 
density - the Local Bubble - initially discovered from the very 
low values of interstellar attenuation of starlight within 100 pc 
of the Sun (Fitzgerald 1968; Lucke 1978; cf. Lallement et al. 
2014a). Lallement et al. (2003) and Welsh et al. (2010) mapped 
the Local Bubble in Na i and Ca ii interstellar absorption, reveal¬ 
ing a central region of very low neutral/singly-ionized atomic 
gas absorption out to a distance of about 80 pc in most Galactic 
directions, and vertically for hundreds of parsecs - presumably 
as far as the Galactic Halo. However, the Ca ii maps did not show 
any sharp increase in absorption at a distance of 80 pc, leaving 
considerable uncertainty about the interface (“wall”) between 
the tenuous local and denser surrounding interstellar medium 
(ISM). It has also been surprisingly difficult to assess the con¬ 
tents of the Local Bubble. Its possible origin in multiple super¬ 
nova explosions (Frisch 1981; Cox & Anderson 1982; Malz- 
Apellaniz 2001; Breitschwerdt et al. 2009) seemed to be cor¬ 


roborated at first by the detection of a diffuse soft-X-ray back¬ 
ground (Snowden et al. 1998). While the putative million-degree 
gas was not seen in the expected extreme-UV emission by the 
CHIPS satellite (Hurwitz, Sasseen & Sirk 2005) nor in O vi 1032 
A absorption (Cox 2005; Barstow et al. 2010), recent works 
making use of soft X-ray shadows, Planck mm-wavelength data 
and Local Bubble maps have now considerably reinforced the 
existence of such hot gas within the cavity (Snowden et al. 
2015). The picture is complicated by the interaction between the 
heliosphere and surrounding ISM, such as charge exchange be¬ 
tween the solar wind ions and interstellar atoms (e.g., Lallement 
et al. 2014b; Zirnstein, Heerikhuisen & McComas 2015; but see 
Galeazzi et al. 2014), indicating local “fluff” of a more neutral 
kind (Frisch 1998). Any novel way or tracer that can be used to 
map the Local Bubble is therefore highly desirable. 

Bring on the Diffuse Interstellar Bands (DIBs). These are 
a set of over 400 broad optical absorption features which are 
ubiquitous in the ISM (Herbig 1995; Sarre 2006). The nature of 
their carriers has been much debated since Mary Lea Heger first 
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observed them almost a century ago (Heger 1922); the current 
consensus is that they are large carbonaceous molecules, with 
some resistance to UV radiation (Tielens 2014). Observations of 
DIBs have shown that they are sensitive to their environment; 
for instance, the ratio of the 5780 and 5797 A DIBs is diagnostic 
of the level of irradiation and/or kinetic temperature (Cami et al. 
1997). Differences in absorption between close binary stars have 
been detected, indicating DIBs trace also tiny-scale structure in 
the ISM (Cordiner et al. 2013). 

In this project we use DIBs to map the Local Bubble, by ob¬ 
serving hundreds of early-type stars within a few hundred pc of 
the Sun. These stars have well determined distances and proper 
motions and future observations could reveal tiny-scale struc¬ 
ture. Also, by probing very local interstellar matter the sight¬ 
lines suffer less from overlapping cloud structures and hence the 
behaviour of the various tracers may be simpler to interpret and 
more telling. DIBs have been used successfully to map extra- 
planar gas (van Loon et al. 2009; Baron et al. 2015) and the 
Magellanic Clouds (van Loon et al. 2013; Bailey et al. 2015) as 
well as parts of the Milky Way (Kos et al. 2014; Zasowski et al. 
2015). The elusiveness of the gas within the Local Bubble, how¬ 
ever, requires very low noise levels on the measurements and 
thus calls for dedicated observations. We have initiated such a 
dedicated observing campaign, first covering the Southern hemi¬ 
sphere presented here, subsequently extending it to cover also 
the Northern hemisphere (Farhang et al. 2015a,b - Papers II 
and III). The combined, all-sky survey is used to apply inversion 
techniques to reconstruct a complete three-dimensional map of 
the density of the DIB carriers (Farhang et al., in prep. - Paper 
IV). In a fifth paper in the series we examine time variations 
of interstellar absorption due to the proper motions of the stars, 
incorporating measurements for stars in common with previous 
works (e.g., Raimond et al. 2012; Puspitarini, Lallement & Chen 
2013). 


2. Measurements - Southern hemisphere survey 

2.1. Targets 

The targets were chosen from the NaiD Local Bubble survey 
of Lallement et al. (2003). The criteria used to select the tar¬ 
gets were that they should be bright stars with B-band magni¬ 
tudes < 8.7, be of spectral type earlier than A5 (to provide a 
clean continuum), have well-known distances from Hipparcos 
and accurately known, high proper motions (for any future time- 
variability studies of the tiny scale structure of the ISM). In all, 
238 stars were observed, with V-band magnitudes ranging from 
2.9 to 8.4. Of these, 109 are of early-A spectral type, 128 of 
spectral type B, and one of spectral type 09.5. Information on 
the observed stars is listed in Table 1, the full version of which is 
available electronically at the Centre de Donnees astronomiques 
de Strasbourg (CDS). 

The locations on the sky are shown in Figure 1. There are 
some regions which have not been sampled as uniformly as oth¬ 
ers; this is partly due to the distribution of potential targets - with 
most of them being concentrated towards the Galactic plane - 
and partly due to adverse observing conditions. To gain a bet¬ 
ter appreciation of the sampling, we show the distributions of 
nearest-neighbour distance in 3D (in pc) in Figure 2; typical 
sampling of 3D space is on 10^0 pc scales. 



Fig. 1. Distribution on the sky of the 238 observed stars as part 
of the Southern hemisphere DIBs survey. The size of the symbol 
indicates the distance from the Sun. 



Fig. 2. Distribution of nearest-neighbour distance, in pc, in the 
Southern hemisphere DIBs survey. 


2.2. Observations 

We obtained long-slit spectroscopic data using the 3.5m 
New Technology Telescope (NTT) at the European Southern 
Observatory, La Silla, Chile. The observations were performed 
in visitor mode on the nights of 2F'-23‘^‘^ March 2011 and 
15*-17* August 2011 (ESO programme 087.C-0405 (A-hB); PI: 
M.Bailey) and 19*-22"‘^ August 2012 (ESO programme 0.89.C- 
0797; PI: M.Bailey). 

The Paint Object Spectrograph and Camera (EPOSC2) were 
used with grism #20 and the 0"3 blue offset slit; this gave a re¬ 
solving power of R = AjAA - 5500, resulting in a two-pixel 
sampled spectral resolution of ss 1.1 A (no binning of pixels 
on the CCD). This combination of grism and offset slit covers 
the wavelength range 5672-6772 A, which includes the major 
45780, 5797, 5850, 6196, 6203, 6269, 6283 and 6614 DIBs. 
Grism #20 introduces a lateral shift of the beam. Therefore, 
the target was placed at column 1680 during acquisition to en- 
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Table 1. Catalogue of observed stars. Only the first three rows are displayed; the full catalogue of 238 objects is made available at 
CDS. The columns list the HD number, equatorial coordinates a and 6 and Galactic coordinates / and b, Hipparcos proper motions 
p and parallax n (van Leeuwen 2007) and distance from the Sun d, V-band magnitude and spectral type (from Simbad - Wenger et 
al. 2000). 


HD 

a (12000) d 

1 

b 

Pa 

Ps 

n 

d 

k 

type 


m ^0 / //^ 

n 

n 

(mas yr ^) 

(mas) 

(pc) 

(mag) 


955 

00 13 55.60 -17 32 43.3 

78.8928 

-77.0855 

17.11 ± 1.05 

4.74 ± 0.45 

3.16 + 0.71 

316 

7.37 

B4V 

2913 

00 32 23.78 +06 57 19.7 

114.5465 

-55.6060 

44.75 ± 0.47 

2.33 + 0.35 

12.35 + 0.55 

81.0 

5.70 

B9.5V 

3580 

00 38 31.85 -20 17 47.7 

98.8792 

-82.5581 

22.24 ± 0.38 

7.68 ± 0.25 

3.97 ± 0.39 

252 

6.71 

B8V 


sure the photons arrived approximately at column 1107 once the 
grism was put in place, which is an area on CCD#40 that is free 
from bad pixels (the dispersion direction runs along the columns 
of the CCD). 

Exposure times were determined with the goal of reaching 
maximum count levels of ~ 40,000-50,000 counts per pixel, 
which is close to but below the CCD saturation level (65,535 
counts); these ranged from 3.4 seconds to 847 seconds. The CCD 
was read out in its fast mode, as the read-out-noise was negligi¬ 
ble compared to the photon noise. A minimum of three expo¬ 
sures were taken for each target to achieve a total signal-to-noise 
level (S/N) of ~ 2,000 after the spectra had been combined. If 
the maximum counts per pixel was < 20,000 a further three ex¬ 
posures were made. Exposure times were altered to accommo¬ 
date changing conditions - the seeing varied between 0!'4 and 
2!'4; some exposures resulted in saturation of the detector dur¬ 
ing moments of excellent seeing as the centering on the narrow 
slit became problematic. 

Before the start of each night a set of calibration frames were 
taken. These comprised a large number of bias frames (> 30), to 
remove the electronic offset level; dome flat fields (~ 60, us¬ 
ing 4x 1,000 W quartz-halogen lamps) to correct for pixel-to- 
pixel sensitivity and throughput variations across the CCD; and 
He-i-Ar arc lamp exposures for wavelength calibration. The large 
number of bias and flat field frames was necessary in order for 
these not to limit the S/N on the calibrated spectra. 


2.3. Data processing 

Standard data reduction steps were applied; creation of master 
frames representing the bias level and flat field, respectively; 
subtraction of the master bias frame from the master flat and 
science frames; division of the science frames by the master flat; 
wavelength calibration on the basis of the He-tAr arc lamp spec¬ 
trum. The spectra were then extracted from the calibrated 2D 
frames using the optimal (weighted) extraction method of Horne 
(1986). This also removed some particularly bright cosmic ray 
impacts upon the detector - though these were few in number 
due to the generally short integration times. 

Two issues were encountered during data processing; bias 
instability and wavelength calibration offsets (for more details 
see Bailey 2014). The fast read-out mode makes use of two am¬ 
plifiers, resulting in different bias levels for both halves of the 
detector. The bias level was found to vary over the course of 
the measurements in an unpredictable manner. It was thus de¬ 
cided in the second and third of the three observing runs to also 
obtain a bias frame after every science observation. The offset 
slit results in part of the CCD remaining unilluminated and we 
therefore decided to use the recorded count level at the edge fur¬ 


thest away from the illuminated part of the CCD, to determine 
an offset value with which to correct the master bias fram^B 

EEOSC2 was originally designed to be mounted at the 
Cassegrain focus (of the ESO 3.6m telescope) where the aper¬ 
ture wheel moved in a horizontal plane inclining only with the 
telescope. Now that EEOSC2 is located at the Nasmyth focus 
(of the NTT), flexure problems are experienced. This is due to 
the fact that the instrument is now turned on its side and gravity 
can influence the position of the slits depending on the rotator 
position of the instrument. The effect of the flexure is a variation 
of the position of the slit on the detector in the dispersion direc¬ 
tion, which results in shifts of the wavelength scale that can be 
as big as several A (EEOSC2 User’s Manual). We corrected for 
these offsets by recalibrating each spectrum using the telluric O 2 
absorption line at T = 6278.1 A. 

2.4. Methodology 

We here first describe the analysis of the 5780 A DIB. A shallow 
broad absorption as well as a stellar absorption line at 5785 A 
were fitted with a second-order polynomial and a Gaussian func¬ 
tion, respectively, at the same time as fitting a Gaussian function 
to the 5780 A DIB. While we considered fitting a Lorentzian in¬ 
stead, the above strategy worked well and is more in line with the 
choice of Gaussian found in the literature (Herbig 1975, Hobbs 
et al. 2008; Jenniskens & Desert 1994; Krelowski & Westerlund 
1988; Westerlund & Krelowski 1988). The wavelength and typ¬ 
ical width of the DIB was taken from the DIB catalogue of 
Jenniskens & Desert (1994), but these were allowed to vary by 
+2 A and 0.5 and 2x, respectively (+2 A and +0.2 A for the stel¬ 
lar line). The instrumental resolution of 1.1 A was sufficiently 
less than the expected width of the 5780 DIB (~ 2 A) to deter¬ 
mine the intrinsic profile shape, but this was not the case for the 
other DIBs under consideration here (~ 0.6-1.2 A) except for 
the 6270 A DIB (~ 3^ A). The actual fitting was coded in idl 
using the mpfit routines, in the range 5773-5788 A. 

The equivalent width EW was determined by integrating the 
Gaussian function; 

EW^J Aexp^- I ^ Ao-V2^, (1) 

where A is the peak intensity, do the central wavelength and cr 
the width (with the Eull Width at Half Maximum, FWHM - 
2.355cr). The errors eo resulting from mpfitexpr were scaled ac¬ 
cording to the deviations from the fit and the degrees of free¬ 
dom / (the number of spectral elements minus the number of 
fitting parameters) as in; 

e = eo (2) 

* A gradient of intensity was observed indicating light still somehow 
made it to the “unilluminated” part of the CCD. 
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Table 2. Overview of the Gaussian fitting constraints, in A. 


Feature 

1 a 




guess 

range 

guess 

range 

5780 DIB 

5780.6 

± 2.00 

1.25 

0.625-2.500 

5797 DIB 

45780+ 16.50 

±0.50 

1.00 

0.500-2.000 

5850 DIB 

45780 + 69.60 

±0.50 

1.00 

0.500-2.000 

He I 

5875.6 

±2.60 

1.30 

0.650-2.600 

Nai D 2 

5889.95 

±2.00 

1.00 

0.500-2.000 

Nai Di 

4d2 + 5.97 

0 

O' D2 

0 

6196 DIB 

45780+415.60 

±0.50 

0.70 

0.350-1.400 

6203 DIB 

45780+422.60 

±0.75 

1.00 

0.500-2.000 

6270 DIB 

45780+489.50 

±0.50 

1.00 

0.500-2.000 

6283 DIB 

45780+503.70 

±0.50 

1.00 

0.500-2.000 

6614 DIB 

45780 + 8 3 3.10 

±0.80 

1.00 

0.500-2.000 


Hence the error in the equivalent width is obtained from: 
e[EW] = ^27r((e[A]cr)2 + (Ae[o-])2) (3) 

The constraints on the fitting to the other absorption lines 
and DIBs are summarised in Table 2. In particular, the wave¬ 
lengths of the other DIBs were guided by that of the strongest, 
5780 A DIB and allowed to vary much less (typically +0.5 A) to 
avoid spurious results. Likewise, the wavelength of the Di (red) 
component of the Na i doublet was tied to that of the stronger D 2 
(blue) component; both components were fitted simultaneously. 
The spectral baseline was fitted with a first-order polynomial 
rather than the second-order polynomial above, except in some 
cases such as the Nat doublet where a shallow, broad absorption 
component was accounted for with a second-order polynomial - 
a second fit was obtained using a first-order polynomial, which 
was divided by the initial fit (including the second-order polyno¬ 
mial and two Gaussian functions) to reveal the broad absorption 
component. The 6196 and 6203 A DIBs were also fitted simul¬ 
taneously and here too a second-order polynomial was required 
for the continuum fitting. A second-order polynomial was also 
required for the 6270 A DIB continuum as well as to take ac¬ 
count of the telluric O 2 line near the 6283 A DIB (before fitting 
the DIB with a Gaussian and first-order polynomial). Finally, 
a second-order polynomial was required to account for mod¬ 
est fringing near the reddest of the DIBs under consideration, 
at 6614 A. 

All spectra and fits were vetted by eye, guided by the value 
for EWjelEW}, and only fits deemed reliable were retained. The 
full table of measurements is available at CDS, with the first 
three entries shown in Table 3. 

3. Results 

3.1. Detected spectral lines and bands 

The He i (D 3 ) line is stellar in origin but the Na i Di & D 2 doublet 
is mostly interstellar. The latter sits on a broader shallow absorp¬ 
tion component, which is a combination of several telluric water 
lines typically one per cent of the continuum level (see Bailey 
2014 for more details). The weak Nai absorption, ~ 5% of the 
continuum on average, illustrates how challenging it is to map 
the ISM within ~ 100 pc around the Sun. 

The DIBs are more challenging still, less than one per cent 
of the continuum. The 5780 A DIB is ~ 0.7% deep on average, 
but clearly detected in many of our targets. The nearby and diag¬ 
nostic comparison 5797 A DIB is < 0.2% below the continuum. 


Table 3. Spectral line measurements. EW stands for equivalent 
width, with e[£’Vk] its uncertainty. Flags / are 1 or 0 for a de¬ 
tection or non-detection, respectively; f - 2 means dubious (ex¬ 
cluded from further analysis). The full, transposed version of this 
table is available on CDS. 


HD 

955 

2913 

3580 

£1^5780 (A) 

0.008 

0.016 

0.020 

e[£lV]5780 (A) 

0.0020 

0.0030 

0.0039 

(£lV/e[£lT])5780 

4.1 

5.4 

5.2 

£1^5797 (A) 

- 

- 

0.008 

e[£lT]5797 (A) 

- 

- 

0.0170 

(£M//e[£lV])5797 

- 

- 

0.5 

£1^5850 (A) 

0.006 

0.011 

- 

e[£W]5850 (A) 

0.0094 

0.0038 

- 

(£lV/e[£lV])5850 

0.7 

2.9 

- 

£lTHe (A) 

0.621 

0.058 

0.015 

e[£lV]He (A) 

0.0180 

0.0061 

0.0104 

(£lV/e[£lT])He 

34.5 

9.6 

1.4 

£lTNaD2 (A) 

0.067 

0.097 

0.022 

e[£'lV]NaD2 (A) 

0.0078 

0.0148 

0.0106 

£lTNaDl (A) 

0.010 

0.043 

- 

e[£lT]NaDl (A) 

0.0046 

0.0109 

- 

(£lT/e[£lT])NaD2 

8.6 

6.5 

2.1 

(£M//e[£lT])NaDi 

2.1 

3.9 

- 

£1^6196 (A) 

0.003 

- 

0.008 

e[£lV]6196 (A) 

0.0021 

- 

0.0026 

£1^6203 (A) 

- 

- 

- 

e[£lV]6203 (A) 

- 

- 

- 

(£iy/e[£lT])6i96 

1.5 

- 

3.1 

(£lV/e[£lV])6203 

- 

- 

- 

£1^6270 (A) 

0.015 

0.002 

- 

e[£lV]6270 (A) 

0.0085 

0.0035 

- 

(£iy/e[£lT])6270 

1.7 

0.5 

- 

£1^6283 (A) 

- 

- 

- 

e[£lT]6283 (A) 

- 

- 

- 

(£lT/e[£lT])6283 

- 

- 

- 

£1^6614 (A) 

0.007 

- 

0.008 

e[£lV]6614 (A) 

0.0066 

- 

0.0044 

(£lV/e[£lT])66i4 

1.1 

- 

1.9 

/5780 

1 

0 

1 

/5797 

0 

0 

0 

/5850 

0 

0 

0 

/ue 

1 

1 

0 

/NaD2 

1 

1 

0 

/NaDl 

0 

2 

0 

/Na dip 

0 

1 

0 

/6196 

0 

0 

0 

/6203 

0 

0 

0 

/6270 

0 

0 

0 

/6283 

0 

0 

0 

/6614 

0 

0 

0 


The 5850 A DIB is very feeble and being next to a much stronger 
stellar feature to the red it is omitted from further analysis. The 
6196 (& 6203 A DIBs are only ~ 0.1% deep on average, but 
sharp and more isolated and thus more useful. The 6270 A DIB 
has a depth of typically a few 0 . 1 % and is safe from the telluric 
O 2 absorption to red, but it is broader and difficult to detect in 
individual targets and hence it is omitted from further analysis. 
Likewise, the 6283 A DIB is omitted; it is but a mere fraction of 
a per cent of continuum and sits on the wing of the strong telluric 
O 2 absorption to the blue. The 6614 A DIB on the other hand is 
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only 0.2% on average but it is well detected often enough to be 
useful for analysis. 


3.2. Correlations 


Figure 3a shows an excellent correlation between the equivalent 
widths of the 5780 and 5797 A DIBs, up to some point beyond 
which increasing strength of the 5797 A DIB is not accompanied 
by a noticeable increase in the strength of the 5780 A DIB. This 
occurs for EW{5191) > 0.04 A, around £1^(5780) ~ 0.15 A. 
Compared to the Friedman et al. (2011) regression fit our sight¬ 
lines seem to probe relatively strong 5780 A DIBs except for 
the sight-lines with the strongest 5797 A DIB. The 5797 A DIB 
favours a neutral environment, whereas the 5780 A DIB favours 
an ionized environment so a strong detection of one should not 
necessarily be accompanied by a strong detection of the other. In 
Table 4 the results are summarised from a linear regression anal¬ 
ysis and a Student t test; the correlation is significant if f > fcnticaU 
where 


t - 


N-2 


(4) 


with correlation coefficient r and N data points {N-2 degrees of 
freedom in this case). These analyses have not taken into account 
the uncertainties of the individual measurements. However, in¬ 
spection of figure 3 shows that either the errors are fairly uniform 
in the cases where the scatter does not vary much along the cor¬ 
relation, or they increase with increasing equivalent widths in 
concert with the scatter. The exact treatment of the uncertainties 
should therefore not have a huge impact as they are already re¬ 
flected in the scatter among the values. Hence we conclude that 
a good correlation is seen when the absorptions are fairly weak 
but the correlation breaks down above EW{5191) > 0.04 A. 

As examples of variations in the EW(5780)IEW(5797) ra¬ 
tio, figure 4 shows the spectra of HD 146284 and HD 169009. 
They have a similar £’VT(5780) but the EW(5797) - 0.058 A of 
HD 169009 is twice the EW{5797) = 0.028 A of HD 146284. 
While HD 146284 is further away from the Sun (d - 315 pc) 
than HD 169009 is = 102 pc), HD 146284 is also further away 
from the Galactic Plane (z - 102 pc) than HD 169009 (z - 2.9 
pc); while the longer column towards HD 146284 through the 
more tenuous extra-planar medium results in relatively strong 
diffuse band absorption, the harsher environment results in a rel¬ 
atively weaker 5797 A DIB (cf. van Loon et al. 2009). 

A good correlation is found between the 6196 and 6614 A 
DIBs (Fig. 3d). This particular pair of DIBs is said to have 
a “near perfect correlation” (McCall et al. 2010). However, at 
the generally low EW values in our study there appears to be 
a change from a tight steep slope up to £1^(6196) < 0.01 A 
and £1T(6614) < 0.02 A, to a shallower relation at greater EW 
values (cf. Table 4). Considering a continuous, monotonic rela¬ 
tion throughout, and the regression fit from McCall et al., there 
is a suggestion in our data that some sight-lines have an anoma¬ 
lously weak 6196 A DIB-e.g., the datum at £1T(6614) = 0.051, 
£W(6196) = 0.006 A (HD 145102). Oddly enough, the McCall 
et al. data also showed some deviations at the weak end but in 
the opposite sense, with the 6614 A DIB being weaker. This ex¬ 
emplifies the bias that can be introduced when sampling differ¬ 
ent sight-lines, as well as distances (longer columns may sample 
multiple clouds hence yielding smoother relationships). 

The 5780 and 5797 A DIBs correlate well with the 6196 and 
6614 A DIB overall (Fig. 3b,c), but the correlations involving 



5770 5780 5790 5800 5810 

wavelength (Angstroem) 



Fig. 4. The 5780 and 5797 A DIBs in the spectra of two stars 
at different distance to the Sun {d) and Galactic Plane (z); Top: 
HD 169009 a.td,z^ 102,2.9 pc; Bottom: HD 146284 a.t d,z ^ 
315,102 pc. The blue spectrum is the original spectrum before 
removing the stellar absorption line around 5785 A; the red and 
green curves are the Gaussian fits to the 5780 and 5797 A DIBs, 
respectively. 


the 5797 A DIB are much tighter than those involving the gen¬ 
erally strong(er) 5780 A DIB (cf. Table 4). This exemplifies the 
different diagnostic power of the 5780 A DIB as compared to 
the other three DIBs - even though we noted that differences be¬ 
tween the 6196 and 6614 A DIB exist too (indeed, there is more 
scatter in the relations between the 5780 and 5797 A DIBs with 
the 6614 A DIB than with the 6196 A DIB). Compared to the 
regression fits from Friedman et al. (2011), the 5780 and 6196 
A DIBs show an identical relation but the 6614 A DIB seems to 
be much stronger in certain (but not all) sight-lines; we confirm 
the Friedman et al. suggestion that the 5780 A DIB is already 
seen when the 6614 A DIB is not (a positive value for a at more 
than 4-cr significance). As an example, figure 5 shows two sight¬ 
lines where the spectra have been superimposed for the regions 
around the 5780, 5797 and 6614 A DIBs; HD 18650 has a larger 
Galactic latitude {b = -25°) than HD 179406 {b = -8°). Clearly, 
their £1T(5780)/£1T(5797) ratio differs by about a factor two, 
but the £VT(5797)/£VT(6614) ratio is very similar. 

The DIBs weakly correlate with the Nat absorption, but 
with a lot of scatter; this indicates that in the typical warm 
neutral/weakly-ionized medium both are seen and they thus each 
correlate with overall gas column density. In figure 6 the be- 
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Fig. 3. Correlations between the 5780, 5797, 6196 and 6614 A DIBs. Distinction is made on the basis of the Na i absorption. Linear 
regression fits are overplotted, including some from the literature. 


Table 4. Linear regression analysis of DIB-DIB correlations of the form Y - a + bX with correlation coefficient r. Where the 
Student t statistic is given, the critical value of t is for a 1% chance that the null hypothesis is true. 


Y/X pair 

selection 


N 

a (A) 


b 

r 

t 

^critical 

verdict 

5780/5797 

all 


46 

0.023 ± 0.005 

2.48 

±0.20 

0.88 

12 

2.7 

significant 

— 

EW{5191) 

<0.04 

39 

0.008 ± 0.005 

3.7 

±0.8 

0.90 

13 

2.7 

significant 

— 

EW(5191) 

>0.04 

7 

0.11 ±0.04 

0.6 

±0.6 

0.44 

1.1 

4.0 

insignificant 

6196/6614 

all 


22 

0.002 ± 0.000 

0.189 ±0.022 

0.89 

8.7 

2.8 

significant 

— 

EW{66U) 

<0.02 

9 

-0.002 ± 0.002 

0.57 

±0.07 

0.95 

8.0 

3.5 

significant 

— 

£IV(6614) 

> 0.04 

13 

0.001 ± 0.001 

0.21 

±0.06 

0.73 

3.5 

3.1 

marginal 

5780/6196 

all 


33 

0.022 ± 0.009 

7.5 

±0.9 

0.82 

8.0 

2.7 

significant 

5780/6614 

all 


32 

0.034 ± 0.008 

1.54 

±0.19 

0.83 

8.2 

2.8 

significant 

5797/6196 

all 


29 

-0.002 ±0.001 

3.0 

±0.3 

0.89 

10 

2.8 

significant 

5797/6614 

all 


25 

0.006 ± 0.002 

0.60 

±0.06 

0.89 

9.4 

2.8 

significant 


haviour of the 5780 and 5797 A DIBs is compared to that of the 
Nat D 2 line. There seems to be a threshold for the 5797 A DIB 
to be seen only above ^^(Nai D 2 ) > 0.1 A; instead, the 5780 A 
DIB can sometimes be seen already at ^^(Nai D 2 ) ~ 0.05 A. 
While this is in keeping with the 5780 A DIB being stronger than 
the 5797 A DIB, the data in figure 6 clearly show the emergence 
of the 5797 A DIB above the threshold whereas the 5780 A DIB 
appears to reach zero strength at zero Nat D 2 absorption. This 
observation may be understood if the 5797 A DIB needs a min¬ 
imum gas column density for the environment to be sufficiently 
cool and/or shielded from UV radiation, whilst the 5780 A DIB 
(and Na i) are seen also in more exposed diffuse ISM. 

Curiously, at moderate ^^(Na i D 2 ) ~ 0.2-0.3 A the DIBs - 
especially the 5780 A one - are particularly strong. Certainly the 


5780 A DIB is not anymore as strong once ^^(Nai D 2 ) ~ 0.4 
A. This could correspond to a greater ubiquity of material within 
which the carriers of the DIBs thrive, at moderate gas densi¬ 
ties, with the 5780 A DIB disappearing at high(er) densities. 
The strong DIB absorption is, however, only seen close to the 
Galactic plane, within about \b\ < 40° (cf. Figs. 6a,b). 

The EW(51%Q)IEW(5191) ratio shows an intriguing be¬ 
haviour with respect to the Nat absorption (Fig. 6c); values 
£’W(5780)/£’W(5797) < 5 are seen across the entire range of 
Nat strength, but values £'VT(5780)/£'W(5797) > 5 are exclu¬ 
sively - and commonly - seen only for ^^(Nai D 2 ) < 0.2 A. 
Smith et al. (2013) find a similar extreme cr-cloud behaviour 
towards ArVelorum (not part of our sample) at a distance of 
~ 165 pc (£W(5780)/£W(5797) > 22 at EW{Nai D 2 ) = 0.12 


6 































































































Bailey et al.: Probing the local diffuse ISM with DIBs. I. 




wavelength (A) 

Fig. 5. The 5780, 5797 and 6614 A DIBs in the spectra of 
HD 186500 and HD 179406; the latter is located closer to the 
Galactic Plane. 

A). This can be understood if the Local Bubble is traced by 
^^^(Nai D 2 ) < 0.2 A, where harsh UV irradiation is preva¬ 
lent. Conversely, we may assume ^^^(Nai D 2 ) = 0.2 A as 
a proxy for the boundary between the Local Bubble and sur¬ 
rounding ISM. The sight-lines with ^^^(Nai D 2 ) < 0.2 A but 
£’VT(5780)/£’1T(5797) < 5 are interesting as these may identify 
compact cloudlets with modest column density but high volume 
density. 

The huge scatter in the above correlation diagrams means 
that Nat is not a very good indicator of DIB carrier abundance 
(Memll & Wilson 1938; Herbig 1995). The reason for this is 
partly due to the ionization potential of sodium; neutral sodium 
will trace weakly ionized gas as well as neutral gas. Therefore, it 
is not a very good discriminator between the two types of envi¬ 
ronment. DIBs, on the other hand, do discriminate between the 
two environments with the 5780 A DIB favouring the weakly 
ionized gas and the 5797 A DIB favouring the neutral gas. The 
use of the NaiD lines is further complicated by saturation ef¬ 
fects already kicking in around ^^(Nai D 2 ) > 0.05 A and the 
Na ID lines not tracing hydrogen column density in a linear fash¬ 
ion (Welty & Hobbs 2001). Hence, DIBs are a more discerning 
probe of the environment and give a more detailed picture of the 
diffuse ISM than Nat can, justifying our aim to map the Local 
Bubble in DIBs (cf. Baron et al. 2015). 

3.3. Maps 

The Galactic latitude distribution of equivalent widths is shown 
in figure 7 for the 5780 and 5797 A DIBs and Na i D 2 . The latter 
shows no relation whatsoever, but the DIBs clearly are stronger 
towards the Galactic Plane, within about \b\ < 40°. The lack 
of dependence of Nat may be due in part to saturation though 
stronger Nat is seen towards more distant objects than in our 
survey (e.g., van Loon et al. 2009). Neither of these three trac¬ 
ers shows clear dependence on distance (colour-coded in the 
graphs). This implies that within a few hundred pc from the Sun, 




D 2 . Distinction is made according to Galactic latitude. 
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• d(pc) > 300 

• 150 < cKpp] S 300 

• 50 < d(pc) S 150 

• d(pc) S 50 





Fig. 7. Galactic latitude distribution vs. equivalent width of the 5780 and 5797 A DIBs and Nai D 2 . Distinction is made according 
to distance. 




Fig. 8 . Distribution on the sky in Galactic coordinates, of the 5780 and 5797 A DIBs and Nat D 2 . The EW values are colour-coded 
and the distance to the Sun is indicated by the size of the symbol (smaller if farther). The directions towards the Small and Large 
Magellanic Clouds and Local Leo Cloud are also indicated. 

small-scale structure dominates the variations between different 
sight-lines. As the DIBs do “detect” the Galactic Plane they seem 
to be more useful tracers to map the Local Bubble than Na i. 

To add a dimension, in figure 8 are plotted the sight-lines in 
Galactic coordinates, i.e. projected on the sky, with EW values 
and distances to the Sun colour-coded and indicated by the size 
of the symbol, respectively. Thus we start to gain a more three- 
dimensional picture of the absorption. In these projections, there 
is a clear gap which corresponds to the Northern hemisphere 
which is covered in Papers II & III. Also indicated are the direc¬ 
tions towards the Magellanic Clouds, which have been studied 
in the past (e.g., Cox et al. 2006; van Loon et al. 2013) and is the 
subject of a large DIB survey presented in Bailey et al. (2015). 

Of more immediate relevance is the location of the Local 
Leo Cloud, a nearby very cold cloud (T ~ 20 K) between 11.3- 
24.3 pc from the Sun. Hence it is well within the cavity of the 
Local Bubble and the closest known cold neutral gas cloud (Peek 
et al. 2011; Meyer et al. 2012). Our targets skirt this cloud but 
some sight-lines from the Northern hemisphere survey crossed 
it. In the present data set two sight-lines just South of the Cloud 
have rather strong 5780 A DIB absorption, without detection of 
the 5797 A DIB. This could be tracing the ionized skin of the 
Cloud. 

Overall, the strongest absorption is concentrated within ~ 

+25° of the Galactic Plane (Fig. 8), but with exceptions. For in¬ 


stance, relatively strong absorption in both the 5780 and 5797 A 

DIBs as well as Nai is seen around I ~ 190° and b -25° to 

-40°. The wall of the Local Bubble may be extending further 
away from the Galactic Plane, possibly if it is relatively close 
to us. On the other hand, sight-lines around I ~ 270° are more 
transparent despite quite long columns - especially when com¬ 
pared to the adjacent sight-lines around / ~ 300°; perhaps the 
Local Bubble extends further in the disc in the former direction. 

The maps in figures 9 and 10 show the EW of the 5780 
and 5797 A DIBs and Nai D 2 , colour-coded, in two projections: 
onto the Galactic Plane and in an orthogonal plane crossing the 
Galactic Centre also (the “meridional” plane). In figure 10 also 
non-detections are shown; in some cases these would mask the 
detections hence the separate plots. The Sun is marked at the co¬ 
ordinates (0,0). The size of the symbols represents the distance 
from the projection plane, with circles representing targets that 
are above/in front and squares representing targets below/behind 
the plane, respectively. 

Areas of strong 5780 and/or 5797 A DIB absorption gen¬ 
erally match areas of strong Nai D 2 absorption and high Na° 
volume density in the Welsh et al. (2010) map. There are fewer 
sight-lines with detected 5797 A DIB absorption, hence these 
maps are sparser; however, including the non-detections, figure 
10 clearly shows the extent of the Local Bubble and in particular 
a “chimney” structure perpendicular to the Galactic Plane, with 
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Fig. 9. Galactic maps of (Left to Right:) 5780 A DIB, 5797 A DIB and Nat D 2 absorption, in projection (Top:) onto the Galactic 
Plane and (Bottom:) in the meridional plane orthogonal to the Galactic Plane. Only detections are included in these maps. 


5797 A DIB absorption confined to the Plane. The 5780 A DIB 
is also seen further out along the chimney, commensurate with it 
being ubiquitious in the extra-planar gas (van Loon et al. 2009). 
So the Local Bubble opens out into the Halo in some tracers but 
not others. Likewise, the sight-lines between I ~ 200° to 270° 
are transparent in the 5797 A DIB but absorption in the 5780 A 
DIB is clearly seen along those same sight-lines. This implies 
translucent density enhancements within the Local Bubble, or 
a fuzzier, extended transition rather than a sharp boundary be¬ 
tween the Local Bubble and surrounding ISM. It is thus impera¬ 
tive not to ignore non-detections in mapping the ISM structure, 
as they clearly trace cavities but also can indicate porosity or 
dumpiness. 

This is emphasized in maps of the ratio of 5780 to 5797 A 
DIB absorption (Fig. 11), where the paler colours trace the neu¬ 
tral gas in the wall of the Local Bubble and the darker colours 
trace the chimney. This is especially clear in the meridional plane 
projection, but even in the Galactic Plane projection the cavity of 
the Local Bubble is clearly characterised by a larger 5780/5797 
A DIB ratio. 

Another exciting thing revealed by this map is that within the 
Local Bubble cavity exist discrete regions of a low 5780/5797 A 
DIB ratio; for instance towards HD 74521 at approx. (-90,-70) 
and (-90 ,h- 70) pc in the Galactic Plane and meridional plane, or 
HD 172910 at approx. (-h 130,0) and (-hl30,-40)pc. These sight¬ 
lines must be hitting clumps of cold neutral gas that somehow 
survive within the harsh environment of the Local Bubble. 


4. Discussion 

4.1. The carriers of the DIBs 

DIBs are currently thought to have molecular carriers, with the 
DIBs representing vibronic transitions from a large number of 
different molecules (McCall et al. 2010). All molecules in a 
given sight-line that give rise to a particular DIB are expected 
to be in their vibronic ground state, implying that two or more 
DIBs arising from that molecule should be part of a vibronic pro¬ 
gression in the molecule; the relative strengths of the vibronic 
bands should then be exactly the same from one sight-line to an¬ 
other (McCall et al. 2010). Therefore, a tight correlation between 
two DIBs is a promising, albeit insufficient, condition to identify 
DIBs with a common carrier (Josafatsson & Snow 1987). The 
problem is, of course, that these molecules may not be in the 
ground vibration state. Differences in excitation of low-lying en¬ 
ergy levels might explain subtle differences between DIBs that 
on the whole behave concurrently. 

The observations presented here of tenuous, generally ex¬ 
posed gas within and around the Local Bubble indicate that the 
6196 and 6614 A DIBs may not form as perfect a pair as recently 
suggested - McCall et al. (2010) determined a correlation coef¬ 
ficient of 0.99, while here we determine r = 0.89 only. While 
our measurements are compatible with those of McCall et al., 
we show clear evidence of differing equivalent width ratios of 
these DIBs and we uncovered a statistically significant break in 
the relation around ^^(bb 14) ~ 0.02-0.04 A. Such behaviour at 
very low column densities may hold clues to the precise relation 
between the carriers of these two DIBs. 
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Fig. 10. As figure 9, with non-detections included. 


The measurements presented here show correlations be¬ 
tween the 5780 and 5797 A DIBs on the one hand, and the 6196 
and 6614 A DIBs on the other. The 5780 A DIB correlates bet¬ 
ter with the 6196 A DIB than with the 6614 A DIB; the 5797 A 
DIB correlates better with the 6614 A DIB than the 5780 A DIB 
does. This suggests that the carrier of the 6614 A DIB is more 
closely related to that of the 5797 A DIB than of the 5780 A DIB. 
The 5797 A DIB is also correlated with the 6196 A DIB, even 
better than the 5780 A DIB is. This suggests that the carrier of 
the 6196 A DIB has similarities to those of the 5780 A DIB and 
- especially - 5797 A DIB. But, as the 5780 and 5797 A DIB 
show different behaviour, the carrier of the 6196 A DIB cannot 
be identical to both the carrier of the 5797 A DIB and that of the 
5780 A DIB. In figure 12a the 6614/6196 A DIB ratio is plot¬ 
ted against the 5780/5797 A DIB ratio. While the 6614/6196 A 
DIB ratio does vary between different sight-lines, overall there 
is no correspondence to the behaviour of the 5780/5797 A DIB 
ratio. The latter is clearly high for low Na° column densities but 
the 6614/6196 A DIB ratio is no different from that in denser 
columns. As in figure 6c, there is also a clear lower limit to the 
5780/5797 A DIB ratio ~ 2, and a similar lower limit is seen for 
the 6614/6196 A DIB ratio. 

Intriguingly, the 6203/6196 A DIB ratio does correlate with 
the 5780/5797 A DIB ratio (Fig. 12b). This is especially con¬ 
vincing for the denser columns, with ^^(Nai D 2 ) > 0.2 A. The 
6203/6196 A DIB ratio seems to drop in more tenuous columns, 
with £W(Nai D 2 ) < 0.2 A (around £1T(5780)/£W(5797) ^ 5); 
perhaps this may be understood if the 6196 A DIB is more read¬ 


ily formed than the 5797 A DIB is (cf. Figs. 3e and 6) - as long 
as the radiation held is not inhibitive. 

We conhrm that the 5780 A DIB favours weakly ionized 
environments and the 5797 A DIB favours neutral environ¬ 
ments, with the relation between the two breaking down above 
EW{5191) > 0.04 A at high statistical signihcance: those denser 
clouds are dominated by neutral matter and the skin no longer 
grows in proportion. Compared to Nat, the 5780 and 5797 A 
DIBs are weaker at higher Galactic latitude. This suggests that 
the carriers are eroded under the particularly harsh extra-planar 
conditions. Also, at \b\ > 40° we could not detect both 6614 and 
6196 A DIBs. 


4.2. Implications for our understanding of the Local Bubble 

The DIB maps of the Local Bubble presented in this work trace 
the cavity and the wall of the Local Bubble rather well. The 
sight-lines with non-detections are of equal importance as the 
sight-lines where DIBs are detected. In this regard the Nai map 
is less discriminative, as Nat absorption is nearly always de¬ 
tected and arises both in neutral and weakly-ionized gas. The 
“DIB-less” sight-lines clearly reveal the extent of the cavity 
within the Galactic Plane, as well as the vertical structure; the 
5797 A DIB map suggests a chimney opening out into the Halo 
though the 5780 A DIB is still detected high up in the extra- 
planar gas. This difference is due to the different alliance of these 
two DIBs to neutral and weakly-ionized gas, respectively. The 
“bursting” of the Local Bubble has been noticed previously (e.g., 
Lallement et al. 2003; Welsh et al. 1999, 2010; Welsh, Sallmen 
& Lallement 2004; Welsh & Shelton 2009; Welsh & Lallement 
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Fig. 11. As figure 9, but for the 5780 to 5797 A DIB ratio. The panels on the right only show ratios derived from detections in both 
DIBs. 


2012) but DIBs may indicate an extra-planar interface between 
the Local Bubble and Halo. 


Within the cavity, taking into account the distances to the 
targets, small neutral cloudlets are found. So even if the Local 
Bubble is filled with hot gas and/or the radiation field is unat¬ 
tenuated, denser, perhaps self-shielded structures do exist within 
it. On the other hand, the “wall” of the Local Bubble is porous, 
with some sight-lines piercing through it out to large distances 
while adjacent sight-lines encounter much absorption. This must 
have implications for the leakage of UV radiation - and possibly 
hot gas - out of the Local Bubble into the surrounding ISM, both 
within the Halo and also the Disc of the Milky Way. 


5. Summary 

In this first paper of a series we have presented the results from a 
Southern hemisphere spectroscopic survey of 238 nearby early- 
type stars, aimed at mapping the ISM within and surrounding 
the Local Bubble in the 5780 and 5797 A DIBs, and Nat. We 
also investigated the relationship with other DIBs detected in 
this spectral region, notably those at 6196 and 6614 A, indicat¬ 
ing subtle differences in their mutual behaviour. The DIB maps, 
including detections as well as non-detections, clearly outline 
the extent of the Local Bubble within the Galactic Plane, and its 
vertical structure opening out into the Halo, and the generally 
harsh conditions within the Local Bubble as indicated by a high 
5780/5797 A DIB ratio throughout. However, small more neutral 
cloudlets were seen close to the Sun as well, whilst the “wall” of 
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Fig. 12. The 6614 to 6196 A DIB ratio vs. 5780 to 5797 A DIB 
ratio. Distinction is made on the basis of the Nat D 2 absorption; 
these four DIBs were never all detected at high Galactic latitudes 
{\b\ > 40°). 


the Local Bubble is far from impermeable. Thus we have demon¬ 
strated the viability of using DIBs to map local ISM, and we will 
use these data in combination with the complementary Northern 
hemisphere survey to reconstruct a full three-dimensional map 
of the Local Bubble. 
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